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Abstract 

We show the existence of self-dual semilocal nontopological vortices in a $^ Chern- 
Simons (C-S) theory. The model of scalar and gauge fields with a SU{2)giobai x U {X)iocai 
symmetry includes both the C-S term and an anomalous magnetic contribution. It is 
demonstrated here, that the vortices are stable or unstable according to whether the 
vector topological mass k is less than or greater than the mass m of the scalar field. 
At the boundary, k, = m, there is a two-parameter family of solutions all saturating 
the self-dual limit. The vortex solutions continuously interpolates between a ring 
shaped structure and a flux tube configuration. 
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1 Introduction 



Topological defects are known to exist in theories with spontaneously broken sym- 
metries. For example, Nielsen and Olesen discovered that Abrikosov type vortices 
(A-N-0) appear as classical solutions of an abelian Higgs model These vortices 
carry magnetic flux but are electrically neutral. Furthermore, for a Higgs po- 
tential, and when the parameters are chosen to make the vector and scalar masses 
equal, minimum energy vortex configurations arise that satisfy first order differential 
equations 0, In this limit, known as the Bogomol'nyi limit, the vortices become 
non- interacting 

In recent years the charged vortex solution |^ of the Abelian Higgs model in (2 + 1) 
dimensions with a Chern-Simons term has attracted a lot of attention in the literature, 
because they can be considered as candidates for anyonlike objects . These vortices 
were shown to exist even in the absence of gauge field kinetic term (Maxwell term) 
[P^. This theory, where the kinetic action for the gauge field is solely the Chern- 
Simons term is known as the pure C-S theory (P-C-S) @, It was recently shown 
that the P-C-S theory with an special choice of the scalar potential V{^) supports 



self-dual topological and nontopological vortices [T^, 0]. For this theory, V{^) is a 
sixth-order potential, and the corresponding Bogomol'nyi limit is obtained when the 
vector and scalar masses are equal. 

More recently an Abelian Chern-Simons model which includes both the C-S term 
and an anomalous magnetic contribution, in addition to the Maxwell term, has been 



studied |]T2[. It was shown that for a special relation between the C-S mass and 
the anomalous magnetic coupling, the equations for the gauge fields reduce from 
second- to first-order differential equations, similar to those of the pure C-S theory. 
Furthermore, it was demonstrated that nontopological charged vortices satisfy a set 
of Bogomol'nyi-type equations for a quadratic potential V{^) = (m^/2)$^, when m 
and the topological masses are equal. This model possess a local U{1) symmetry, so 
we will refer to it as the local $^ model. 

On the other hand, new interest in the study of string-like defects has arisen af- 



ter the observation made by Vachasparti and Achucarro [|13| that the Nielsen- Olesen 
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vortex solution can be embedded into a larger theory which has a global SU{N) sym- 
metry in addition to the local U{1) symmetry, these objects are known as semilocal 



vortices. Semilocal vortices also appear in Chern-Simons theories; Khare WM ob- 
tained semilocal vortex solutions when the local $^ Abelian Higgs model of Refs. 
[|T0| , [TT[] is extended to a semilocal one. 



In this paper we consider a SU{2)g ® U{l)i model, where only the overall U{1) 
phase is gauged and SU (2) is a global symmetry. We will consider a simple $^ scalar 
potential, so we are necessarily in the symmetric phase of the theory; nevertheless we 
shall find static minimum energy nonotopological vortex configurations. We will be 
interested in the stability of these semilocal vortices. We will be able to write a lower 
bound for the energy (Bogomol'nyi bound), when the scalar (m) and the topological 
(k) masses are equal. The lower bound is saturated when the fields satisfy a set of 
first order differential equations (self-dual or Bogomol'nyi equations). The self-dual 
vortices are neutrally stable, but it will be shown that the vortex solutions become 
stable or unstable according to whether k is less than or greater than m. The semilocal 
$^ model exhibits a richer spectrum of solutions as compared with the local model. 
In particular the vortices can assume both a ring shaped structure (typical of the 
local P-C-S models) as well as a flux tube structure ( typical of the A-N-0 vortices). 
We should also point out that the stability condition m > k is exactly opposite to 
the one obtained for the semilocal Nielsen-Olesen vortices 



The paper is organized as follows. In Sec. 2 we present the semilocal model and 
discuss their general properties. In Sec. 3 we derive the self-duality (Bogomol'nyi) 
equations. The derivation of this self-duality equations is much more involved than 
the original Bogomol'nyi derivation for the Nielsen-Olesen model, therefore we present 
it with some detail. Sec. 4 is devoted to the study of the cylindrically symmetric 
vortex solutions. First, it is shown that the local $^ vortex solution can be embedded 
into the larger semilocal string; then a generalization of the ansatz is presented, this 
generalized solution represents a two-parameter family of solutions all saturating the 
self-dual limit. Section 5 deals with stability analysis of the solutions. Sec. 6 presents 
a numerical study of the self-dual solitons. Concluding remarks comprise the final 
section. 
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2 Model 



The model consists of a complex doublet 

*=ftM. (2.i; 



^2 . 

with only the overall phase gauged. The semilocal Lagrangian is 

£ = -\f^,F^^' + '^e^^^A^F,^ + - , (2.2) 

where F^y = d^Ai, — d^A^, we use natural units % = c = 1 and the Minkowski-space 
metric is g^j^^ = diag{+l, —1, —1); fi = (0, 1, 2). The covariant derivative includes 
both the usual minimal coupling plus the anomalous magnetic contribution: 

D^^ = {d^ - leA^ - i^e^y^F^'^)^ , (2.3) 



with g the anomalous magnetic moment |T^. Notice that the Lagrangian Eq. ( p.2|) 
has a global SU{2) symmetry and a local U{1) symmetry. We should also remark 
that it is an specific feature of a (2 + 1) dimensional world, that a Pauli-type coupling 
[i.e., a magnetic coupling) can be incorporated into the covariant derivative, even for 
spinless particles |12[. In fact, it was demonstrated in Ref. that radiative 



corrections can induce a magnetic coupling for anyons, that is proportional to the 
fractional spin. 

The equations of motion for the lagrangian in Eq. (|2.2|) are 



Df,D^'<^ = -m^<!> , (2.4) 



e^^uad'^lF'^ + ^ J"] = Ju- kF, . (2.5) 

The last equation has been written in terms of the dual field, = \e^yaF^°', and 
the conserved matter current = [p, J) is given by 



■^M 2 



(2.6) 
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The energy momentum tensor is obtained by varying the curved-space form of the 
action with respect to the metric 



1,„ ^,9 rn?,^,') 

- Va$ $ 

2I A I 2 I I 



+ H.c. , (2.7) 



where = — ieA^ includes only the gauge potential contribution. Notice, that 
both the Chern-Simons and linear terms in g do not appear explicitly in T^^,. This is 
a consequence of the fact that these terms do not make use of the space-time metric 
tensor g^i,; consequently, when g^^ is varied to produce T^u no contributions arise 
from these terms 0. 

For the selected scalar potential the theory is always in a symmetric phase. In 
this case the theory possess two propagating modes in the trivial sector (excitations 
around the vacuum): an scalar field excitation with mass m and a massive vector 
mode with mass k. Instead, if we would consider a broken symmetry phase the 
theory would contain three propagating modes, because the gauge field acquire two 
distinct masses due to the P and T violating terms |T3, IT^ . 



3 The self-duality equations 

There is a particular relation between the C-S mass and the anomalous magnetic 
moment for which the Eq. ( ^.5| ) for the gauge fields reduce from second- to first-order 
differential equations |]TB|, O, ED], similar to those of the P-C-S type M. To get this 



limit notice that if the following relation holds 

2e , , 

K = , 3.1 

9 

then it is clear that the Eq. (|2.5|) are solved identically if we choose the first order 
ansatz 

F, = -J,, (3.2) 
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that have the same structure as the equations of the P-C-S theory 0. We will refer 
to the previous conditions as the P-C-S limit. However, we should notice that the 
explicit expression for differs from the usual expression of the P-C-S theory, because 
according to Eq. (p.6|) and Eq. (|2.3|) receives contributions from the anomalous 
magnetic moment. This P-C-S equations (Eq. ( p.2|) ) imply that any object carrying 
magnetic flux must also carry electric charge (Q), with the two quantities related 



as Q = — In what follows we shall work in the limit in which Eq. (|3.2|) and 

Eq. (|3.1| ) are valid, so we consider Eq. (|3^) as the equation of motion for the gauge 
fields, instead of the Eq. ( |2.5| ). 

In the so called Bogolmol'nyi limit all the equations of motion are known to 
become first order differential equations 0; furthermore, it is possible to write the 
equations of motion as self-duality equations. In reference [0, it was shown that 
for the local model, it is possible to find Bogomol'nyi-type equations for a quadratic 
potential (m^/2)$^, when m and the topological masses are equal. In that case, a 
rotationally invariant vortex ansatz was substituted in the expression for the energy 
functional, then the Bogomol'nyi equations were found in terms of the reduced number 
of functions that appear in the ansatz. However, the self-duality equations in terms 
of the original fields $ and without assuming the rotational invar iance were no 
presented there. In order to discuss the vortex solutions in the extended semilocal 
model we require to know the form of the self-duality equations. The derivation of 
this self-duality equations is not as straightforward as in the original Bogomol'nyi 
derivation for the Nielsen-Olesen solutions, therefore we present with some detail the 
deduction of the self-duality equations for the present model. 

In order to obtain the self-duality equations of motion let us consider the energy 
density (Too) in the static limit written in terms of the magnetic B = ^eijF^^ and 
electric Ei = fields, this yields: 



Tf 



00 



1 - 



B^ + \E\^] + 



ivo$r+(Vi$)*Vi$ 



-I- -m |$| 



(3.3) 



We recall that the covariant derivative = — ieA^ includes only the gauge 
potential contribution. Now, we notice that we can exploit the /i = component of 
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the equation Eq. ( p.6| ) together with the C-S equations of motion Eq. ( p.2| ) to express 
^0 in terms of the magnetic field 



An 



B 



From this equation we can observe, that for a time-independent vortex solution, the 
Aq component can not be set to zero otherwise the magnetic flux would vanish. Using 
the above equation, the first and the third terms in Eq. (|3.3| ) can be added together 
to obtain the simplifying result 



1 
2 



1 



1 



(3.5) 



Likewise, from the fi = i components of Eq. (|3.2| ) and Eq. (|2.6| ) we can express 
the electric field in terms of V,$: 



1 



K 



(3.6) 



In order to simplify the second and fourth terms of the energy density Eq. ( p.3| 
it is convenient to define a new covariant derivative Di according to the relation 



L) . = Vi + iHFi = di- ieAi + iHFi , 
where we have introduced the auxiliary function H defined as 



(3.7) 



H 



K 



e|$| 



21^12 



1 



em 



e2|,^l2 



'1 



(3.8) 



Using the previous results (Eq. (|3.6|) and Eq. (|3.71) ), we can now find that the second 
and fourth terms in Eq. ( |3.3|) add together to give the result 



1 - 



l$l 



|E|2 + (v,<l')*V,$ = (A$) A$, 



(3.9) 



where summation over latin indices is from i = 1 to i = 2. 

Employing the results of Eq. ( |3.5| ) and Eq. ( |3.9| ) we can write down the energy 
E = J Tood'^x as 
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^ = (s^iif (i - $1*1') + \ (A*)' A* + 5"'i*r) ■ (3.10) 

The energy written in this form is similar to the expression that appears in the Nielsen- 
Olesen model. Thus, starting from Eq. ( |3.1CI| ) we can follow the usual Bogomol'nyi- 
type arguments in order to obtain the self-dual limit. The energy may then be 
rewritten, after an integration by parts, as 



E 



d X 





(1- 


e2 






e2 




2 



21^12 



1 - 



1/2 



+ 1 [DiTiDi]^ 



2 



± — i dl-n ± -<p dl-A, 



6 J r=oo 

where the vectors Q and A are defined according to the relations: 



(3.11) 



A,; 



1 - 



1/2 



A, 



K 

e2|<|)|2 



1 - 



1/2N 



(3.12) 



For any nontopological soliton the asymptotic conditions are such that $ — at 
spatial infinity. Thus, the line integral of A in Eq. ( p.ll| ) vanishes, whereas the line 
integral of Q yields the magnetic flux 



dl-Q 



dl-A=<^ 



B 



(3.13) 



In what follows, we consider only those configurations that fulfill the condition |$| < 
n/e. Consequently, from Eq. ( |3.11| ) we can conclude that the energy is bounded 

2 

below; for a fixed value of the magnetic flux, the lower bound is given hy E > 
provided that the potential is chosen as a with the critical value m = i.e. 

when the scalar and the the topological masses are equal. Therefore, in this limit we 



8 



are necessarily in the symmetric phase of the theory. From Eq. ( p.ll| ) we see that the 
lower bound for the energy 



E = -\^b\ = -IQI, (3.14) 

e e 

is saturated when the following self-duality equations are satisfied: 

bi^ = izDa'^', (3.15) 



B = ±- ^T^, (3.16) 



1 



e2|«i»|: 



where the upper (lower) sign corresponds to positive (negative) value of the magnetic 
flux. We should remark that these self-duality are valid both for the local as well 
that for the semilocal $^ model. The second of these equations implies that the 
magnetic field vanish whenever $ does. For the local $^ model |jT2[ the finiteness 
energy condition forces the scalar field to vanish both at the center of the vortex 
and also at spatial infinity; consequently for the local model the magnetic flux of the 
vortices lies in a ring, rather than being concentrated at the center as in the A-N-0 
solution. As we shall see below, for the present semilocal model $ does not vanish at 
the origin; furthermore, the vortices can assume both a ring shape as well as a flux 
tube form. 

With regard to the self-duality equations (Eq. ( |3.15| ) and Eq. ( p.l6| )), we note 



that they are similar to those found in other models, but there are some important 
differences. First, we have to point out that in order to derive the self-dual limit it was 
essential to introduce a new covariant derivative (Eq. p.7| )) that is different from 
the original (Eq. ( |2.3| )) that appears in the Lagrangian. Of particular interest is to 
compare the present results with those of the Nielsen- Olesen (A-N-0) and the pure 
Chern-Simons (P-C-S) models, {i) Similar self-duality equations arise both in the 
local and semilocal |l^, |22| A-N-0 models for a potential of the form (|$P — v'^)'^ 



when the scalar and vector masses are equal. For the self-dual A-N-0 model the 
Eq. (^) is similar, but with the normal covariant derivative = — ieA^ instead 
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of D^] furthermore, Eq. ( p.l6| ) is replaced by one of the form 5 ~ — and the 
magnetic field is maximum at the center of the vortex, {ii) Self-dual soliton solutions 
have also been found in a local IITUI O and semilocal versions of the P-C-S with 



no Maxwell term and no magnetic moment contribution for a sixth order potential 
of the form — f^)^, when the parameters are chosen to make the vector 

and scalar masses equal. Again the form of Eq. ( |3.15|) holds, but with the normal 



covariant derivative = — ieA^] whereas, Eq. ( |3.16|) is now replaced by one of 
the form B ~ — for this kind of vortices the magnetic flux is localized 

within a ring around the center of the vortex. 

4 Vortex solutions 

The local $^ vortex solution of reference I1T2| can be embedded in one of the com- 



ponents of $. Explicitly for the semilocal model a static rotationally symmetric 
configuration of vorticity n can be written in in terms of the following ansatz: 



A{p) = -9^ , Aop =-% 

ep e 



Hp) = ^/(P)exp(-m0) I I , (4.1) 

where (p, 6) are the polar coordinates. If the previous ansatz is substituted into the 
self-duality equations of motion (Eq. (|3.15| ) and Eq. ( |3.16|) ) the equations reduce ex- 




actly to those of the local vortices [|I^ , and hence we have self-dual nontopological 
solitons also in the present semilocal model. That this is the case, may be seen by 
noting that if the lower component of $ is set to zero the resulting model is exactly 
the local $^ model. This means, that a solution of the local $^ model is automatically 
a solution of the semilocal model. Most of the properties of the semilocal string are 
identical to those of the local vortex: namely they represent charged flux tubes that 
carry fractional spin as well as magnetization. However, there may be a difference: 
the stability of the semilocal vortex may be different from the stability of the local 
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one. Indeed, if one holds the lower component of $ to zero and perturbs only the 
upper component, it is exactly the same as perturbing the local $^ solution that 
is known to be stable when m > k, However, we can also perturb the lower 



component of $, so we should check whether the semilocal vortex remains stable for 
m > K. 

Before proceeding to study the stability of the semilocal vortices we notice, follow- 
ing Hindmarsh that the ansatz Eq. (|4.1| ) can be generalized. In fact, the most 



general ansatz which maintains the cylindrical symmetry is 



A{p) = -e ^^^^ , Ao(p) = -%), 



K 

e 



(4.2) 



/ f{p)exp{—in9 
\F{p) exp {—imS 

The same as in reference it is sufficient to examine only the case m = 0; this is 
equivalent to add a cylindrically symmetrical perturbation to the lower component of 
the ansatz Eq. ( [4.1[ ). With this ansatz the self-duality equations of motion Eq. ( |3.15D 
and Eq. ( p.l6|) become 



r h ' 

J r (/2 + F2) V h 

(a-n)F npF / 1\ , , 

where we have introduced the dimensionless variable r = /tp, primes denote differen- 
tiation with respect to r and we make use of Eq. ( p.4| ) and Eq. ( p.l6| ) to solve for the 
p = component of the gauge field that yields 



h=[l-f'-F'Y. (4.4) 
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The boundary conditions are selected in such a way that the solution Eq. ( ^.21 ) is 
non-singular at the origin and gives rise to a finite energy solution; then, the problem 
is to solve Eq. ( [4. 3D subject to the following boundary conditions: 



a = n, / = 0, F' = 0, at r = , 

a —a, / ^ 0, F — >• 0, as r — > oo . (4.5) 

Notice that at spatial infinity the value a(oo) = —a is not constrained. Consequently, 
we shall see that for nontopological solitons the magnetic fiux is not quantized, but 
rather it is a continuous parameter describing the solution. 

The system of differential equations Eq. (|4.3| ) looks rather involved, however the 
system can be simplified considerably if we notice that the equations for / and F are 
not independent. In fact it can be easily demonstrated, that subject to the boundary 
conditions Eq. (|4.5|), / and F are related according to the relation 



± 



/, (4.6) 



where tq is an arbitrary parameter. Exploiting this result, the three differential 
equations (Eq. ( [4.3|) ) reduce to two; it is convenient to write them in terms of the 
functions a(r) and h{r), this brings the equations to the form 



a' = r 



h' 



jh' - 1) 
h 

jh' - 1) 



n 



1 + 



2n 



(4.7) 



In what follows, we select the signs (upper signs in al the previous equations) corre- 
sponding to positive magnetic flux {n > 0). The equations for n < are obtained 
with the replacement a — > —a, f f, F F and h — > —h. This previous system of 
equations is more amenable to be treated numerically and various properties of the 
solutions can be discerned by general considerations. In particular the large distance 
behavior of the solutions yields 
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-a + ^—r—, + 0(r-^"+^) 

(a - 1) r2"-2 ^ ' 



h = l-^ + 0(r-^"), (4i 



where C is a constant. From these asymptotic expressions we see that the magnetic 
field, B = a'/r, falls off like r~^". Likewise the electric field \E\ oc h' falls off like 
^-{2a+i)^ We should remark that although the theory includes gauge fields with mass 
K, the magnetic field departs from the usual e~'^^ asymptotic behavior and becomes 
a power law. From the asymptotic expression for a(r) we notice that there should 
be lower bound on the values of a. Indeed, a > 1 so the second term in the first of 
Eq. ( [4.8|) is subleading compared with the first. For small r a power- series solution 
gives 



n 



2ho 
(hl- 



2hl 



2ho 



(4.9) 



The constant is not determined by the behavior of the fields near the origin, but 
it is instead a parameter of the vortex solutions. The value of ho may be restricted 
by requiring the proper behavior (Eq. ([4.5|) ) as r ^ oo. 

Once that the boundary conditions are known, the topological numbers of the 
soliton can be explicitly computed. With the ansatz Eq. (|4.2| ) the magnetic field is 
B = a'/r, and using the boundary conditions Eq. (|4.5|) the magnetic fiux yields 



f 271 2tx 

= cPxB = — (a(0) - a(oo)) = — (n + a) 
J e e 



(4.10) 



271 . , 2tt 
— (a(0) - a(oo)) = — 
e e 

notice that for nontopological solitons the magnetic fiux is not quantized. The solu- 
tions are also characterized by the charge Q, spin 5* (which is general fractional) and 
magnetic moment M: 
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y = —K(pB = n + a , b = —{a — n 

e ^ 

M = / r^—dr. 4.11 

e Jo dr 

Notice that the magnetic flux, the charge and the spin can be calculated explicitly, be- 
cause they depend only on the boundary conditions. Whereas, the magnetic moment 
depends on the structure factor of the vortex configuration. 

5 Vortex stability 

In the self-dual limit the soliton energy is proportional to its charge (Eq. ( |3.14]) ). 
This result is somehow similar the one obtained for the Q-balls although in our 
model the soliton solutions are time-independent, whereas the Q-balls are necesarily 
time-dependent. For nontopological solutions the soliton charge is of the same type 
(Noether-charge) as that carried by the elementary excitations of the theory . Conse- 
quently, we should check wether the soliton is stable against emission of elementary 
particles. As we shall see at the self-dual point (m = k) the vortices are neutrally 
stable. For other values of the parameters the semilocal vortices are stable when 
m > K, and are unstable otherwise. Notice that this situation is reversed as compared 
to the results obtained for the semilocal Nielsen- Olesen (A-N-0) vortices. Indeed, 
Hindmarsh demonstrated that the A-N-0 vortices are stable when the mass of 



the scalar particle (m) is smaller that the mass of the vector field (m„). We shall 
explain the reason of this difference at the end of the section. 

As stated above, the self-dual vortices are neutrally stable. This fact can be 
easily demonstrated on account of the relation between the energy and the charge 
(Eq. ( |3.14| )): E = K,\Q\/e. First, we recall that the mass of the elementary excita- 
tions of the theory (scalar particles) is m and the charge e. Because of the charge 
conservation a decaying soliton should radiate Q/e "quantas" of the scalar particles. 
Therefore, the energy of the elementary excitations is £^ = mQ/e. This indicates 
that the vortices are at the threshold of stability against decay to the elementary 
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excitations, due to the fact that the ratio En/£ = n/m is equal to one at the critical 
point m = K. 

In order to discuss the soliton stability we need to discuss the existence of two 
sums rules that can be derived from Eq. ( [4.7|) . To obtain these sum rules we follow 
the reasoning introduced by Khare |2^. The first sum rule can be established by 



multiplying the first of Eq. ( [4.7|) by r, integrating over dr and using the asymptotic 
values of a, this yields 

nco r (I - h'^) 

a + n= -dr. (5.1) 

Jo h 

For the second sum rule we multiply the first Eq. (|4.7] ) by a, after integrating and 
using the second of Eq. (^^ ) we obtain 

a^-n^ = 1n\ r^—dr + 2 (l-h'^)rdr. (5.2) 

io i^(^r_y''dr Jo ^ ^ ^ ' 



From these two sum rules we an derive an upper limit for a. On account of Eq. ( [4.4| 
we see that < h < 1; thus we can combine Eq. (|5.1|) and Eq. (|5.2|) to obtain 



r (1 - h^) , r°° / , 9\ , 

a + n = / ^ — ; -dr > r (1 - h ] dr , 

Jo h Jo ^ ^ 

From Eq. ( [4.7D we notice that for the positive flux solution da/dr is negative; then, 
the previous equation yields the upper limit a < n + 2. As mentioned in the previous 
section, a is also bounded below: a > 1. Consequently, we find the bounds 1 < a < 
n + 2. This result should be compared with the bounds that were found for the local 
$^ model: n < a < n + 2 PT| . In the limit rg ^ Eq. (|5.2|) implies the lower limit 



n < a, that is the same bound that was found found for the local vortex |^ . 

We now take up the discussion of the the vortex stability for solution away from 
the self-dual point {m ^ k). We shall follow |^ and and consider the effect of 



perturbing the potential around the self-dual point. We write the potential as 
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ni$i) = ^^^i$r, (5.4) 

that corresponds to a scalar mass m = k, + e, with |e| ^ 1. The functions h, /, F 
and a will be modified from their self-dual configurations. We write them as: 



a{r) = asd{r) + eai{r) + 0{e' 

f{r) = fsd{r) + eMr) + 0{e'), 

F{r) = F,rf(r) + eFi(r) + 0(e2), 

h{r) = hsd{r) + eh{r) + 0{e^) . (5.5) 

From now on, we use the subindex sd to denote the values of the quantities in the self- 
dual point. Therefore, the functions asd{r) and hsdij") are solutions of the differential 
equations Eq. (|4.7| ) and they satisfy the sum rules Eq. (|5.1| ) and Eq. (|5.2|) . On 



substituting these expressions into the ansatz Eq. ( |4.2|) the resulting expressions for 



the energy functional Eq. ( p.ll| ), to 0(e), can be recast in the form 



E^ = -\Q\ + e^ r(l-hl)dr. 5.6 

Notice that the charge (Eq. ([4.11|) ) is modified as compared to its self-dual value, be- 
cause the asymptotic value a changes a = + where ai is unknown. Neverthe- 
less, when one consider the ratio of the soliton energy to the energy of the elementary 
excitations {£ = niQ/e) we obtain that to order e, ai does not contribute: 



-i = l-e—- {n + asd)- r (l - hi) dr . (5.7) 
£ eQsd I Jo ^ ' \ 

We notice that on account of the sum rule Eq. ( |5.1|) the quantity inside the square 
brackets is positive (see Eq. ( ^TsD ). Thus, we conclude that when e > (m > k) the 
soliton solution is stable against dissociation into free scalar particles. Whereas, the 
soliton is unstable for e < (m < k). 

As stated before, the stability condition is exactly opposite to that obtained for 
the semilocal A-N-0 vortices @]. The reason is traced to the fact that, although the 
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semilocal A-N-0 vortices are not strictly topological, they inherit some topological 
properties of the local model; e. g the boundary conditions are not modified if we 
vary the parameters of the model and the magnetic fiux is quantized. Instead, in 
the model both the local and semilocal vortices are of nontopological nature. 
Using a variational procedure Hindmarsh found that the A-N-0 vortex energy 
increases above the Bogomol'nyi limit when m > rriv] but stability requires that the 
Bogomol'nyi limit must be saturated a condition that contradicts the possibility of 
having finite value for the magnetic fiux.. In the case of the semilocal <I>^ vortices 
Eq. ( ^.6|) seems to suggest that the soliton energy also increases above the Bogomol'nyi 
limit when m > but this is not necessarily true, because now the charge is not 
quantized and therefore its value is also modified. Furthermore, for nontopological 
solitons the relevant quantity is the ratio En/S of the vortex energy to the energy of 
the elementary excitations. For values of the parameters such that m > n the mass of 
the elementary excitation increases and the net effect is that the ratio En/S decreases 
(Eq. (p. 71)), so it renders the soliton stable. 

6 Numerical solutions 

To complete the analysis of the previous sections we present the numerical solutions of 
the self-duality equations (Eq. {\i.7f) ) with the boundary conditions given in Eq. ( [4.5|) . 



The solutions depend on two parameters ro and ho (see Eq. ( ^Iq) and Eq. (|4.9| )) 
that define the regular solutions. Regular solutions exists in certain region of the 
parameter space [hovsro); first we recall that < ho < 1. Furthermore, the values 
of ro and ho should be selected in such a way that the boundary conditions Eq. ( [4.5|) 
at r = oo can be met. In this regard our results may be summarized in Fig. 1, where 
we show that there are three regions in parameter space. In region (///) there are no 
acceptable solutions consistent with the boundary conditions at infinity (Eq. ( [4.8| )). 
Both in regions (J) and (//) there are acceptable vortex solutions, but the properties 
of these solutions vary from one region to the other. In region (/) the vortices have 
a fiux tube structure with the magnetic field peaked at the origin; whereas in region 
(//) the vortices have a ring structure, in this case the maximum of the magnetic field 
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is not at the center of the vortex. We explain these results below, in what follows we 
consider the case of isolated vortex i.e. n = 1. 

The origin of the boundary between regions (J) and (//) can be easily understood 
from the expression for the asymptotic behavior of the fields at small r Eq. (|4.9| ). The 
field h{r) starts at the value ho at the origin and should always approach the value 
h 1 at spatial infinity. According to Eq. ([4.9|) the slope of h{r) at the origin can 
either positive or negative. Looking at Eq. ( [4.9|) it is convenient to define the critical 
value he as 

h. = ^ r ■ (6.1) 

The boundary between regions (/) and (//) corresponds to the equation ho = hdro). 
In region (J) we have h^ < he, the slope of h{r) at the origin is positive and the 
function increases monotonically from ho at r = to the value h = 1 at spatial infinity. 
Instead, in region {II) h^ > he the function h{r) decreases from /iq until it reaches 
a minimum value hmin, and from this point it increases to the value h{oo) 1. 
It is clear that this kind of solution will be acceptable if hmin > 0, otherwise the 
equations become singular, see Eq. ( [4.7[ ). The boundary between regions (//) and 
(/J J) represents the points at which /i^m = 0; the values at this boundary are found 
numerically. In Fig. 2 we plot h{r) for a selected tq = 1 and several values of the 
parameter /iq. For the selected value of tq = 1 and looking at Fig. 1, we expect to 
have acceptable solutions for any value of /iq; however, the behavior of the solutions 
should change as ho cross the point h^ = 0.4142. For ho < he the solution increases 
monotonically, whereas for ho > he the solution develops a minimum away from the 
origin. These properties of the solutions are verified in the plots presented in Figs. 1 
and 2. The corresponding solutions for the function a(r) are presented in Fig. 3. As 
expected from Eq. ( [4.7| ) the function a{r) decreases monotonically from a(0) = 1 to 
the value a = —a at r = oo. 

In Fig. 4 we show profiles for the magnetic field B = a'/r as a function of r for 
several values of ho and ro = 1. It is useful to notice that using Eq. ( |3.16| ), Eq. ( [4.2| ) 
and Eq. ([4.4|) the magnetic field can be written in terms of the field h{r) as 
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e n[r) 

Thus, as ho increases the magnetic field at the origin (5 oc (1 — ho) /ho) decreases. 
Furthermore, as explained above, in region (J) we have ho < h^ and the field h{r) 
increases monotonically with r; consequently the magnetic field decreases monotoni- 
cally and the vortex has a fiux tube structure. However as the value of ho increases 
above he the magnetic field smoothly diminishes at the origin and develops a maxi- 
mum away from the origin. Therefore in region (//) the magnetic field of the vortex is 
localized within a ring around the normal core. We recall that from the Chern-Simon 
equation (Eq. ( p.2| )) the charge distribution is proportional to the magnetic field. 

In Fig. 5 we show the electric field {E oc h') for the values of the parameters 
mentioned above. We notice that in region (/) the electric field has a maximum at 
finite r ( the minimum is reached at the boundaries) , however for parameters in region 
(//) the field develops also a minimum away from the origin. 

The vortex solutions are characterized by the quantum numbers that were calcu- 
lated in section 4, namely magnetic fiux (Eq. (|4.1CI| )), charge and spin (Eq. (|4.11|) ). 
These quantities as well as the soliton energy (Eq. ( p.l4| )) depend on the asymptotic 
value of the field a(oo) —a. The value of a is in general a function of ro and ho- 
We recall from section 5 that a is bounded 1 < a < 3 ( for n = 1). Fig. 6 shows the 
behavior of a as a function of ho for a fixed value of tq = 1; we find that a approaches 
the upper limit a — 3 as /iq — ^ 1- Instead, in Fig. 7 we present the plot of a as a 
function of tq for ho = 0.5; in this case as ro ^ oo we notice that a approaches the 
lower limit a ^ 1. For rg < 0.69 there is no vortex solution, in agreement with the 
results of Fig 1. 
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7 Conclusions 



By extending the local model of reference to a SU{2) f/(l) semilocal model 
we found time independent charged vortex solutions. It was first necessary to derive 
the self-duality equations of motion (Eq. ( p.l5D ) in terms of the original scalar ($) 
and vector fields (Af^). This self-dual limit occurs for a simple $^ potential, when 
the scalar and the vector topological masses are equal. Although it is not possible to 
find exact solutions of the self-duality equations, many of the elementary properties 
of the vortices can be uncovered. The self-dual vortices are neutrally stable, but 
including small perturbations away from the self-dual point it was demonstrated that 
the vortices are stable or unstable according to whether the vector topological mass k, 
is less than or greater than the mass m of the scalar field. It is interesting to remark 
that this stability condition is oposite to the one obtained for the semilocal A-N-0 
vortices [^. 

The properties of these vortices are indeed remarkable. Both the local and 



semilocal [0, ^ A-N-0 vortices are known to have a flux tube structure with the 
magnetic field peaked at the center of the vortex. For the charged vortices of the local 
P-C-S models the magnetic field is concentrated within a ring surrounding the center 



of the vortex with the magnetic field vanishing at this point [|T0|, |TT], |T9|. Instead 
in the semilocal versions of the $^ P-C-S model |jT^ the magnetic field is nonzero 
at the center of the vortex. In the present semilocal model the vortex solution can 
have both a ring shaped structure or a flux tube structure. The solution continuously 
interpolates between the two configuration as we move from region (J) to region (//) 
in parameter space (Fig. 1). 

There are several points to explore. In particular the study of the complete de- 
scription of the multisoliton solutions and the behavior away from the self dual limit 
deserve special attention. Of particular interest would be to study the properties of 
Chern-Simons vortices upon quantization, because they can be considered as candi- 
dates for anyonlike objects in quasiplanar systems. 
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FIGURE CAPTIONS 

FIG 1. The crossover lines separating the three regions in parameter space for the 
n = 1 vortex solution. The boundary between region (J) and (//) is obtained from 
ho = /ic(^o) where is given in Eq. ( |6.1D . In region (///) there are no acceptable 
soliton solution, because hmm becomes negative. 

FIG 2. A plot of h{r) as a function of r for a fixed tq = 1 and values of the 
parameter Hq of 0.25,0.40,0.50,0.75,0.95. For a vortex of vorticity n = 1. 

FIG 3. A plot of a(r) as a function of r for a fixed tq = 1 and values of the 
parameter Hq of 0.40,0.75,0.95. For a vortex of vorticity n = 1. 

FIG 4. The magnetic field in units of K^/e for the n = 1 vortex solution, a fixed 
ro = 1 and values of the parameter ho of 0.25,0.40,0.50,0.75,0.95. 

FIG 5. The electric field in units of n/e for the n = 1 vortex solution, a fixed 
tq = 1 and values of the parameter ho of 0.25,0.40,0.50,0.75,0.95. 

FIG 6. Behavior of a as function of ho for ro = 1 and n = 1. 

FIG 7. Behavior of a as function of tq for ho = 0.5 and n = 1. 
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